Irregular orbits generate higher harmonics 
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The spectrum of higher harmonics in atoms calculated with a uniformized semiclassical propagator 
is presented and it is shown that higher harmonic generation is an interference phenomenon which 
can be described semiclassically. This can be concluded from the good agreement with the quantum 
spectrum. Moreover, the formation of a plateau in the spectrum is specifically due to the interference 
of irregular, time delayed, trajectories with regular orbits without a time-delay. This is proven by the 
absence of the plateau in an artificial semiclassical spectrum generated from a sample of trajectories 
from which the irregular trajectories (only a few percent) have been discarded. 

PACS numbers: 32.80Wr, 5.45Mt, 3.65Sq 
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The generation of higher harmonics (HHG) is an in- 
triguing and experimentally well confirmed phenomenon 
which results from the non-linear response of a micro- 
scopic system to a strong laser field pLE| . HHG has been 
studied in simple but illustrative models numerically and 
analytically ||-|(|, for reviews see . Thereby, two 
striking features have been identified, namely the occur- 
rence of a "plateau", i.e. the almost constant intensity 
of the harmonics over a wide range of orders N, and the 
sharp "cutoff" at a certain maximum order N max of har- 
monics. These features have been explained in terms of 
a simple quasiclassical argument 

A closer inspection, however, reveals that only the cut- 
off can be explained with this argument that involves a 
phase matching condition for the semiclassical amplitude 
imposing constraints on the actions of representative clas- 
sical orbits. In the case of an initially bound electron one 
obtains the intuitively appealing picture that the electron 
must return to the nucleus in a certain time correlated 
with the period (frequency) of the laser field to generate 
higher harmonics |Q . If the electron has too much energy 
(which it would need to generate extremely high harmon- 
ics) it is too fast to fulfill the matching condition. Hence, 
the matching condition does explain the cutoff, or more 
precisely, it predicts that the conditions for HHG are un- 
favorable for N > N max . On the other hand this does 
not explain the existence and origin of the plateau for 
A" < N max since the cut off condition does not provide 
a reason why the probability for HHG should be (almost 
uniformly) high for A~ < N max as it is found in exper- 
iments and in numerical simulations. Indeed, only in 
quantum simulations is the plateau found, classical sim- 
ulations do not yield a plateau. This raises the question 
whether the plateau is due to inherently quantum me- 
chanical effects, such as diffraction or tunneling, or if it 



is a pure interference phenomenon that can be explained 
semiclassically. 

In order to answer this question one must carry out a 
full semiclassical calculation of HHG which has not be 
done so far. This is probably due to considerable techni- 
cal difficulties since the chaotic dynamics of the explicitly 
time dependent problem renders a standard semiclassical 
treatment (even for one spatial degree of freedom) in the 
framework of the van Vlcck propagator || impossible. 
However, using a uniformized propagator following the 
ideas of Hermann and Kluk (l(],[ll| we have succeeded 
in obtaining a converged semiclassical spectrum of HHG. 
Moreover, we are able (i) to prove that HHG is a pure 
interference effect, and (ii) to identify the different types 
of trajectories which interfere with each other. 

We have performed our calculation with the "canon- 
ical" model system for the interaction of a strong laser 
field with a one-electron atom, described by the Hamil- 
tonian 



H =p 2 /2 + V(x) + E x cos ut 



(1) 



where V(x) 



-(x 2 



- 1 / 2 with , 



2 is the so called 



"soft core" potential (atomic units are used if not stated 
otherwise). With this choice of a the ground state en- 
ergy in the potential V corresponds to that of hydrogen, 
E = —1/2 a.u.. The other parameters which will be used 
are E = 0.1 a.u. and u> = 0.0378 a.u.. We propagate a 
wavepacket ^(x,t) according to 



\9{x,t)) = U{t)\*(x,0)). 



(2) 



The initial wavepacket has its center xq — Eq/loq = 70 
atomic units away from the nucleus (located at x = 0) 
and is defined as 



*(sb,0) = ( 

7T 



1/4 



exp 



2 XXQ 



(3) 
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FIG. 1. Quantum (a), classical (b), and semiclassical (c) 
spectrum of higher harmonics according to Eq. (Q) . 

with A a b = a — b, 7 = 0.2236 a.u. and po = a.u.. Under 
these conditions of a scattering experiment the cutoff for 
HHG occurs at 2U P = Eq/2lu 2 , see also || where the 
same initial conditions have been used apart from the 
width 7 which does not occur there. 

The Gaussian form of ij)(x, 0) allows one to express the 
semiclassically propagated wavefunction in closed form 
as an integral over phase space [Et| , 



= Jp^fiS J j d( l d P R y(Pt,qt) exp (^S{p t ,q t ) 



exp 
exp 



2 xqt + n PtA 



(4) 
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where S(pt, qt) is the classical action of a trajectory at t, 
and 



R-f(Pt,qt) 



1 



M qq + M pp -ij 2 ?iM, 



1/2 



(5) 



is composed of all four blocks M a ^ — d 2 S / {dadb) of the 
monodromy matrix. 

From the time dependent wavefunction we construct 
the dipole acceleration 
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FIG. 2. Examples for direct (dashed line) , stranded (dotted 
line) and trapped (solid line) trajectories, see text. 



dV(x) 



d(t) = - (#(t) 
from which the harmonic spectrum 
a{w) = ( d{t)e iut dt 



*(*) 



& (6) 



(7) 



is obtained by Fourier transform. Typically 10 6 trajec- 
tories are necessary to converge d(t) from Eq. (||). For 
comparison we have also determined d(t) quantum me- 
chanically (Fig. |j) using standard Fast Fourier Transform 
split operator methods (FFT) to compute the wavefunc- 
tion ^>(x. t). 

Figure [TJ demonstrates that a plateau and a cutoff 
are visible in the quantum (a), and in the semiclassi- 
cal (c) harmonic spectrum, but not in the classical (b) 
one. Since the semiclassical spectrum (b) and the quan- 
tum spectrum (a) are very similar we may conclude that 
HHG can be described semiclassically. Furthermore, the 
absence of the plateau in the classical spectrum (b) sug- 
gests that it is due to an interference effect of different 
types of classical trajectories contributing to the semi- 
classical result (c). 

Among the classical trajectories from which the semi- 
classical dipole acceleration Eq. (g) is constructed we can 
distinguish trajectories which suffer a time delay when 
passing the nucleus (i.e. x ~ 0) from the "mainstream" 
trajectories which are not slowed down. Furthermore, 
among the time-delayed trajectories we can identify two 
groups. 

Trajectories of the first group (dotted line in Fig.|2|) get 
"stranded" on top of the barrier of the effective potential 
V e ff(x) = V{x) — EqX. The second group is formed by 
trajectories which become temporally "trapped", (solid 
line in Fig. |^). The trapped trajectories are chaotic in 
the sense of an extreme sensitivity to a change in initial 
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FIG. 3. Deflection function qt(pi) for t = 3T and Xi = xo = 70 a.u. demonstrating the chaotic character of the trapped 
trajectories. The arrows indicate the range of the next higher enlargement. 



conditions. This is clearly seen in the deflection function 
Itipi) (Fig. H) where the final position q t of a trajectory 
at fixed time t is plotted versus its initial momentum pi 

& 

One sees that in certain intervals of pi small changes 
in pi lead to a completely different qt with the result that 
the deflection function exhibits a fractal structure. The 
fractal initial conditions (for a fixed final qt) belong to 
those trajectories which are trapped in the potential for 
a certain dwell time (the solid lines in Fig. ||). 

The time-delayed irregular orbits are responsible for 
the higher harmonics since their contributions interfere 
with those from the mainstream trajectories. The in- 
terference manifests itself in a dephasing in the dipole 
response d(t) of Eq. (Q) after the first encounter with the 
nucleus (roughly after the time t = T = 2tt/lo for our ini- 
tial conditions) as can be seen on Fig. ||b. At this time 
the peak at about Pi ~ — Ah a.u. , emerges in the deflec- 
tion function, see Fig. [| This corresponds to the return 
of the nucleus in the case of an initially bound electron 
as discussed, e.g., in The rich structure of this peak 

emerges for longer times (see Fig. ||) necessary to resolve 
the fractal dynamics on a fine scale of the initial condi- 
tions pi. The dephasing in d(t) is clearly an interference 
phenomenon since it does not occur in the classical dipole 
response (Fig. ||a). 

Having identified the orbits, or equivalently, the ini- 
tial conditions, which are responsible for the higher har- 
monics we can artificially construct a harmonic spectrum 
without those contributions to double check that they are 
really responsible for HHG. This has been done in the 
semiclassical spectrum of Fig. ||b where the time-delayed 
trajectories (about 3% of all initial conditions) have been 
discarded. Clearly, the plateau has disappeared and the 
spectrum is similar to the purely classical spectrum with 
trajectories for all initial conditions included (Fig. |l|b). 



Discarding only the trapped trajectories (0.6%) smears 
out the cutoff and leaves a reduced plateau for lower har- 
monics (Fig. ||a). Hence, the quantitative semiclassical 
reproduction of the quantum HHG spectrum together 
with the absence of higher harmonics in the classical case 
(Fig. |l|) and in the semiclassical case if irregular, time- 
delayed trajectories are discarded (Fig. ||) confirms our 
explanation of the origin of the higher harmonics. 
To summarize we have shown that higher harmonic 
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FIG. 4. Classical (a) and semiclassical (b) dipole accelera- 
tion according to Eq. (0) . 
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FIG. 5. Semiclassical HHG spectrum as in Fig. |l| but with- 
out trapped trajectories (a), and without time-delayed trajec- 
tories (b), see text. 

generation can be interpreted as a semiclassical interfer- 
ence effect between regular and time-delayed trajectories 
of the electron. The time-delay is either due to a tempo- 
ral trapping which generates chaotic dynamics or due to 
a stranding on top of the potential barrier. Along with 
this time delay goes a characteristic difference in action 
compared to the undelayed mainstream orbits. Analyt- 
ical quasiclassical approximations of various kinds have 
been used to derive this phase difference which can ex- 
plain the cutoff ^j^. However, as demonstrated here, 
the full semiclassical expression is far more complicated 
since for the HHG spectrum chaotic trajectories exhibit- 
ing a fractal deflection function are essential. The chaotic 
character of the irregular orbits allows them to have a rel- 
ative large effect in comparison to their weight among all 
initial conditions (of the order of 1%) because their in- 
stability leads to a dramatic increase of their weight i? 7 
in Eq. (p|) in the course of time. This increase makes an 



accurate semiclassical computation rather difficult. Re- 
markably, despite the chaotic dynamics of the trapped 
trajectories, one can obtain a converged semiclassical 
spectrum if a proper semiclassical propagator such as the 
Hermann-Kluk propagator is used which does not break 
down at the (abundantly occuring) caustics. The result- 
ing semiclassical harmonic spectrum agrees well with the 
quantum spectrum. 
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